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Abstract 

In high T c superconductors a wide ranging connection between the doping dependence of the 
transition temperature T c and the room temperature thermopower Q has been observed. A "uni- 
versal correlation" between these two quantities exists with the thermopower vanishing at optimum 
doping as noted by OCTHH (Obertelli, Cooper, Tallon, Honma and Hor). In this work we provide 
an interpretation of this OCTHH universality in terms of a possible underlying quantum critical 
point (QCP) at T c . Central to our viewpoint is the recently noted Kelvin formula relating the ther- 
mopower to the density derivative of the entropy. Perspective on this formula is gained through 
a model calculation of the various Kubo formulas in an exactly solved 1-dimensional model with 
various limiting procedures of wave vector and frequency. 
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I. INTRODUCTION 



Universal properties of strongly correlated matter are particularly interesting, since 
they hold the promise of revealing the fundamental physics of these systems. An exam- 
ple in the case of heavy fermion systems is the well-known Kadowaki-Woods relation [l| 
A 7~ 2 ~ 1.0 x 10~ 5 ufi cm between the specific heat coefficient 7 and A, the T 2 coefficient 
of the resistivity. Understanding the origin of this universal value has led to considerable 
theoretical progress [2J. In the case of High T c materials, Obertelli, Cooper and Tallon 



(OCT) [3| observed that the thermopower for several cuprates vanishes in the vicinity of 
optimal doping. Honma and Hor (HH) [4J extended their analysis, which was based on the 
phenomenological scale 1 — = 82.6(x — 0.16) 2 [5|, by first showing that the thermopower 
for all the cuprates collapses onto a universal curve as a function of the doping level with 
a zero crossing at a hole content of x ~ 0.23. Since this universal scale is valid for all the 
cuprates, they then advocated [4] that it can be used as an independent calibration of the 
doping level. Consequently, the mapped out the phase diagram of the cuprates in which 
x was determined from the thermopower not the Presland [5] scale. Using this scale, they 



showed [4] that optimal doping for 19 of the 23 cuprates studied corresponded with the zero 
crossing of the thermopower. The exceptions are 4 single-layer materials. Assuming this 
is not a coincidence, it is reasonable to conclude that the mechanism of high-T c and the 
vanishing of the thermopower js] share a common origin. 

The observation of OCTHH has stimulated considerable thought in the community [7 10]. 
In the context of the Hubbard model, some js, 9] have argued that dynamical spectral weight 
transfer ceases at the doping at which the thermopower vanishes, thereby defining a quantum 
phase transition (QPT) where the upper-Hubbard band decouples and Fermi-liquid theory 



obtains. While cluster calculations on the Hubbard model [ll|] support this interpretation, 
no simple model has been studied in which exact statements can be made regarding the 
vanishing of the thermopower and the onset of a quantum phase transition. The present 
work is stimulated by this situation, and we present below two key ideas and a set of model 
calculations that provide a natural framework for understanding such a universality. 

The first key idea in our work is the interpretation of the thermopower as being largely 
determined by thermodynamics, rather than transport aspects, such as velocities and relax- 
ation times [2 1 1 . While this idea of thermodynamic domination of thermopower cannot be an 
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exact statement, it leads to the Kelvin formula proposed by Shastry and coworkers [lo[ Q 



in the spirit of Lord Kelvin's original treatment 



13(. The Kelvin formula for thermopower 



Qk is obtained by computing the slow limit of an exact formula at finite q,u, and is given 



by 



22 1 



q e \dxJ TV 

where S is the entropy density, x is the density of carriers in the system and q e is the charge 
of the carriers (— |g e | for electrons) |23| . In brief, this approximate formula for thermopower 
captures the enhancements due to all fluctuations that influence the thermodynamics of 



a many body system. While theoretical benchmarks of this approximation exist [10], it 
is also useful to check its consequences directly for the high-T c systems. Using standard 
thermodynamics, one gets the following relation between the temperature dependence of 
the Kelvin thermopower and the specific heat variation with particle density x: 

9Qk 1 / dC v \ (dl 



dT T \ dx J vt t^o \dx y ^ , 



(2) 



where 7 is the low-temperature coefficient of the specific heat. This equality comprises 
a relationship between two independent experiments and can therefore be tested with ex- 



perimental data. Fig. [T] plots the left and right side of (2) based on thermopower (3] and 
electronic specific heat data [ji], 16] for Bi 2 Sr 2 CaCu 2 0s+5 and Tl 2 Ba 2 Cu06+«5. While the 
agreement is far from perfect and deviations are certainly expected since the correspondence 
between the left and right-hand sides of Eq. (T5]) is only expected to hold at T = 0, the signs 
and orders of magnitude of each side of ([2]) are compliant with the data, suggesting the level 
of accuracy one can expect from the Kelvin formula in these systems. 

The second key idea involves an underlying quantum critical point jsl, (QCP) as the 
origin of the OCTHH universality. We suggest that optimum doping corresponds to a QCP 
at T = 0, so that at room temperature, the entropy has a maximum as a function of the 
density. From the Kelvin formula, this feature would explain why the room temperature 
thermopower changes sign at optimal doping. 

These ideas may be illustrated within a simple fermionic model exhibiting a QCP. The 
reader is forwarned that this model is quite unphysical because the number of particles is 
not conserved; nevertheless, its exact solvability makes it invaluable for the purpose at hand. 
It also has the essential feature that the parameter tuned to reach the QCP is thermody- 
namically conjugate to the average particle number density. Because the fingerprints of the 
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FIG. 1: A test of the reliability of the Kelvin formula, as explained in the text, for Tl2Ba2Cu06+5 
at T = 100 K (left) and BSSCO at T = 120 K (right). The red solid curves are dq/dx and the 
blue dotted curves are \q e \ dQ/dT. These curves would coincide in the T — > limit if the Kelvin 
formula were exact . 



T = transition at finite temperature include a maximum in the entropy as a function of 
the density, the properties of the model lead naturally to the OCTHH universality 



II. MODEL 



We analyse the anisotropic quantum XY model in the presence of a transverse field in 
one dimension, which, after a Jordan-Wigner transformation, is described by the fermionic 
Hamiltonian Q'Qj 

7 = - J2 [ c i c '+! + rc H+i + h - c ] - h X^ 1 - 2c i°J- ( 3 ) 

i i 

Here — 2h is the dimensionless chemical potential. Due to the T terms, the number of 
fermions is not conserved and hence only the average particle number may be fixed. The 
Hamiltionian can be diagonalized via the Bogoliubov transformation: its eigenvalues are 

e k = ±2y/(h - cos k) 2 + T 2 sin k 2 , - vr < k < vr. (4) 

For any value of T, the system has a QCP at h = 1, where the spectral gap A = 2|(1 — h)\ 
vanishes continuously. The T-h phase diagram is shown in Figure 2a. There are two 
low-T phases. For h < 1, the low-T phase has an energy gap with a high density of 
fermions. At T = 0, T > and with \h\ < 1, the equal-time spin correlation function 



C(R,t = 0) = (of of , fl ) [25| in the limit i? — >■ oo is non-zero, being equal to [18] C(oo,0) = 
2r^n [r 2 (l ~~ h 2 )] 1 ■ This result indicates the presence of magnetic long-range order in the 
ground state. For h > 1, the low-T phase has a spectral gap and with a low number density. 
It corresponds to a quantum paramagnet in the original spin model. 

The lines A = t shown in Fig. [2h have vanishing excitation energies and represent the 
crossover between the low-T gapped phases and the intermediate phase, where the dimen- 
sionless temperature is defined as t — k B T/J. The average particle density x = (qCj) 
is 

1 dk tankC-Ar 21 ) • (5) 



1 

X = 2 



7T 



It ' \e{k)\ 

The critical particle density x c (T) = x(h — 1, T) at T — is plotted in the inset of Fig. [2^. 
Due to the nontrivial mapping between h and x in (151), the t—x phase diagram is dependent 
on the magnitude of T, as is shown in Fig. [2b. On the other hand, the t-h phase diagram is 
independent of T. The lines corresponding to a vanishing gap at h — 1 bend as one moves 
along the t axis. The bending angle varies with T and is larger for smaller values of T. In 
addition, the crossovers A = t between the low-T and high-T phases move as T is varied. 

The entropy S for this model may also be found exactly 17|. For all t, S is a maximum 
as a function of density at x(h — 1), as shown in Fig. [3j The location of the maximum has 
its origins in the large number of micro-states possible for the original spin system at h = 1. 
The locus of the maxima of S depends upon V and is depicted as the central line in Fig. [2b. 
We term this the "peak line" . The peak line has not received attention in earlier works but 
plays an important role for our purpose. 



A. Choice of current operators 

Eq. ([3]) is the spin-less version of the ID BCS reduced Hamiltonian. As the total number 
of particles is not conserved, the standard continuity equation for the charge density p(i) 
does not hold and there is some ambiguity as to the choice of local charge current operator. 
We work with the standard charge current operator given by J{n) = i(c4c n+ i — c\ l+l c n ) 
(here q e = 1). With this choice, the continuity equation for the charge density has pair 



sources and sinks as in the BCS problem [19j. The continuity equation for the local energy 
density is standard 9H 9 ^ + (Je (n+i) — Je (n)) — with energy current operator given by 
Je (n) = i (1-r 2 ) (cl l _ 1 c n+1 - c^ +1 c n _i ) +i 2h fc+Cn-i - cj^Cn) +i 2hT (c^c], - c„c n _i J . 
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FIG. 2: a) Phase diagram of the ID quantum model described by the Hamiltonian ([3]) as a function 
of the chemical potential h and temperature t. There is a quantum phase transition at t = and 
h = 1 for any value of anisotropy parameter T. The entropy is a maximum along the central vertical 
line that we call the "peak line". The two lines A = t mark the crossovers from the low-T gapped 
phases to a high-T phase. Inset: The particle density at the critical point x c = x(h = 1,T = 0), 
which varies significantly as a function of the anisotropy T. b) t— (x — x c ) phase diagram of the 
same model for various values of the anisotropy parameter F. The peak line and the crossover lines 
in t—h plane split into many lines in the t—x plane, depending on the value of V. 



ca 



CO 



0.14 
0.12 
0.1 
0.08 
0.06 
0.04 
0.02 




[a] T=0.25 

k B T=0.05J 



=0.1 J 
=0.1 5J 
=0.2J 




0.05 0.15 0.25 0.35 0.45 
X 



CD 



if) 



[b] 



r=i 



0.05 
0.04 
0.03 
0.02 
0.01 




0.05 0.15 0.25 0.35 0.45 
X 



FIG. 3: Plots of the entropy S as a function of hole density x at several values of t for a) T = 0.25 
and b) r = 1. Notice the clear shift in the location of the maximum of S as t increases for T = 0.25. 
This corresponds to the "peak line" shown in Fig. [2b. 



It is interesting that both the charge and the energy currents are conserved in this integrable 
model. However, the calculations do not simplify on account of this feature and so we will 
not pursue its consequences further. 
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B. Formulae for the thermopower 



For a generic ID system there are four possible linear response formulae for the 
momentum- and frequency-dependent thermopower: 

rft(^) - (6a) 

TQ 2 (q,uj) = -— (6b) 

TQ 3 {q,uj) = (6c) 

Xp( q ),JHq){U) 

T Q 4 (q,uj) = — 6d 

Xj(,),jt(,)(w) 

Here A(q) = Y2 n e~ iq ' n A(n) denotes the Fourier transform of the local operator A(n). It can 
be shown by integration by parts that these four formulas are equivalent if [p(0), H(0)] = 0, 
i.e., particle number is conserved. The fast limit of (I6bj) yields the celebrated Kubo formula; 
the slow limit of (16a]) is ([T|) the Kelvin formula jl^j ]. 

The thermopower was calculated in the slow and fast limits for all four linear response 
expressions. The Kubo formula is given by 

n _ n fast _ l E fc sin(fc)^F(fc,^) 

QKub ° Q2 T J2 k sin(k)Mk)^ ' (7j 

where F(k,h) = 2[(h - cos k)A x {k) + T sm^B^k)], A x {k) = 2 %°; (fc)) , B 1 (k) = 2 -^^, 
and Uk is the Fermi function. Fig. |4] plots the thermopower calculated from the Kubo and 
the Kelvin formula for t = 0.1 and T = 0.75. Note that the Kubo formula result changes 
sign across the QCP through a divergence. It is positive for x < x c and negative for x > x c . 
The Kelvin formula captures the broad features of the thermopower from the Kubo formula 
result, but it instead changes sign through a zero. 

The fast limit of Q\ and the slow limit of Q2 are 0(g) and hence vanish. Aside from 
those deriving from Qi and Q2, there are four more possible expressions for the thermopower 
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FIG. 4: The thermopower, as calculated from the Kubo formula (blue dotted curve) and the Kelvin 
formula (red solid curve). In both the cases, the thermopower is positive for x < x c and negative 
for x > x c . The result shown is for t = 0.1 and T = 0.75. 

which come from the slow and the fast limit of Q% and Q^. 
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(8) 



Here V(k) = 2sin(fc), V H (k) = 2hV(k) - 2(1 - T 2 ) sin(2fc) and e'(k) = 2fM. Fig. [5] shows 
plots of the formulae (JE]). These four expressions for the thermopower give very similar 
results: they are positive for all values of the particle density x and do not show any sign 
change across the QCP. 

In the end, one must make a decision as to which thermopower is definitive. On the basis 
of linear response theory, one knows that formulae derived from fl6a|) and (16bj) concern the 
response of the system to a longitudinal electromagnetic field while formulas ( J6cl) and ( I6dl) 
concern a transverse electromagnetic field. Since the thermopower measured experimentally 
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FIG. 5: The thermopower, as calculated from the formulae Both the slow and the fast limits of 
Q3,4, which give the thermopower in the presence of a transverse electromagnetic field, are positive 
for all values of the particle density x and show no sign change across the QCP. The plot shown is 
for t = 0.1 and T = 0.75. 

is actually the response to a longitudinal electromagnetic field, we believe that only results 
from ( l6a| and ( 16bj) contain information about the physical thermopower j^ . 

III. CONCLUSIONS 



Supposing that the choice of thermopower is correct, the OCTHH universality is partially 
vindicated by the model. The Kelvin and Kubo formulae for the thermopower exhibit a 
sign change at x c (T) « x c (T = 0), in parallel to the case in the cuprate superconductors. 
Inherently, the model studied here is much too simple to be taken seriously as a microscopic 
model of a real material: in particular, it can do nothing to correlate optimal doping with a 
QCP. It does, however, illustrate how a sign change in the thermopower might ultimately be 
connected to a QCP. It also illustrates how an equilibrium construction such as the Kelvin 
formula may be used to approximate the behavior of a transport quantity such as the Kubo 
formula. Finally, we point out that the method used by Vidyadhiraja, et al. llj to locate the 
quantum critical point relies on the maximum in the entropy and hence is closely linked with 
the Kelvin formula. The fact that their state-of-the-art calculations pinpoint optimal doping 
with a maximum in the entropy represents an independent corroboration (albeit not an exact 
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statement) that the sign-change in the thermopower (that is the OCTHH universality) does 
signify a QCP which has been further connected to the Mottness collapse jsl, I9], Q. 
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